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Abstract 

Acoustic  signals  of  small  arm’s  fire,  the  muzzle 
blast  and  the  shock  wave  generated  by  a  supersonic 
bullet  in  air,  are  difficult  to  mask  and  can  be 
exploited  for  localization  of  the  hidden  sniper.  The 
paper  presents  the  system  of  acoustic  measurements 
based  on  a  number  of  both  directional  and 
omnidirectional  microphones  detecting  the  shock 
wave  only,  yielding  exact  solution  for  the  sniper 
direction  in  spite  of  certain  measurement  errors  in 
the  directional  measurements.  The  system  has 
a  self-correcting  ability  concerning  the  sound 
directional  measurements  which  contributes  to  the 
system  technical  feasibility.  Auxiliary  muzzle  blast 
measurements  would  yield  the  sniper  position. 

1.  Introduction 

Acoustic  sniper  localization  exploits  either  of  two 
acoustic  phenomena:  (1)  the  muzzle  blast  detected 
by  a  system  of  microphones,  which  indicates  the 
rifle  position  [1],  and  (2)  the  shock  wave  generated 
by  a  supersonic  bullet  which  produces  a  very 
characteristic  acoustic  signal  [2]  detected  by  wide¬ 
band  microphones.  Two  directional  measurements 
of  the  shock  wave  cone  (the  sound  propagation 
direction  and  its  arrival  time  at  a  given  microphone 
position)  suffice  for  evaluation  of  the  bullet  path, 
assuming  that  its  velocity  is  constant.  The  time 
delay  of  the  muzzle  blast  signal  helps  to  pinpoint 
the  sniper  position  [3]. 

Only  the  acoustic  signals:  the  muzzle  blast  to  some 
degree,  and  primarily  the  shock  wave  [4]  generated 
by  a  supersonic  bullet  passing  by,  could  be 
exploited  for  the  effective  sniper  detection  as  they 
cannot  be  in  principle  masked.  It  requires  the 
appropriate  system  to  be  installed  in  the  protected 
area  that  makes  necessary  acoustic  measurements 
and  evaluates  the  sniper  position  or,  at  least,  his 
direction,  allowing  one  to  direct  a  proper  counter¬ 
fire. 


Such  systems  are  studied  in  recent  literature  [3,  5], 
where  the  two  above  mentioned  acoustic 
measurements  are  used  to  localize  the  sniper. 
Assuming  the  known  bullet  velocity  and  correct 
detection  of  the  muzzle  blast  among  possible  nearly 
simultaneous  other  blasts,  these  two  measurements 
allow  one  to  compute  the  sniper  position  with  rather 
low  accuracy;  experiments  show  the  50%  chance  of 
successful  counter- fire  [6].  The  substantial 
difficulty  arises  from  at  least  two  problems 
connected  with  the  acoustic  measurement  of  the 
muzzle  blast:  1)  it  is  the  relatively  low-frequency 
signal  whose  arrival  time  can  be  detected  with  low 
resolution,  and  2)  it  propagates  over  rather  large 
distance  in  air  of  generally  variable  properties 
bending  the  sound  propagation  path  and  somehow 
enlarging  the  propagation  time.  It  contributes  to  the 
evaluation  error  of  the  sniper  distance. 

The  interesting  shock  wave  signal  propagates  over 
a  rather  small  distance  trough  assumed  reasonably 
uniform  air.  Moreover,  it  is  the  very  characteristic 
N-shaped  signal  with  the  rising  time  much  below 
lps  thus  the  shock  wave  arrival  time  can  be 
detected  with  high  resolution  by  a  proper  wide  band 
microphone  measuring  the  high-frequency  signals. 
It  makes  the  acoustic  measurements  of  the  shock 
wave  much  more  accurate  and  reliable  than  the 
measurement  of  the  muzzle  blast.  Even  the 
information  about  the  bullet  dimension  -  thus  on  the 
used  sniper’s  arm  -  is  included  in  the  shock-wave 
N-shaped  signal,  and  can  be  exploited  for  defense 
purposes. 

These  peculiar  properties  of  the  shock-wave 
measurements,  whose  significance  is  seemingly 
under-  estimated  in  the  existing  literature,  deserve 
high  attention  in  the  considered  problem  of  the 
acoustic  sniper  localization.  In  fact,  only  the  shock 
wave  arrival  time  measurements  are  discussed  in 
the  paper.  Naturally,  using  the  shock  wave 
measurements  alone,  one  cannot  evaluate  the 
distance  covered  by  the  speeding  bullet  and  the 
position  of  the  sniper.  The  muzzle  blast 
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measurements  must  be  used  for  that  according  to 
the  well  developed  method  presented  in  [5]. 
However,  the  shock-wave  measurements  are  able  to 
yield,  as  it  is  shown  in  the  paper,  the  accurate 
direction  of  the  bullet  velocity,  assuming  its 
closeness  to  the  sniper  direction,  that  is  sufficient  to 
direct  the  counter- fire.  Moreover,  the  accurate 
direction  to  sniper  allows  us  to  detect  the  sniper 
with  other  than  acoustic  observations,  the  optical 
detection  of  the  hot  muzzle  gases,  for  instance. 

The  evaluation  of  the  bullet  velocity  direction 
requires  measurements  of  the  shock  wave  arrival 
times  at  a  number  of  observation  points 
(microphones).  As  it  is  known,  the  shock  wave 
generated  by  a  supersonic  bullet  has  the  form  of 
a  cone;  its  axis  is  the  bullet  path  (assumed  straight) 
and  the  conical  angle  that  depends  on  the  bullet 
velocity  (assumed  constant). 

The  measurement  data  must  be  sufficient  to  solve 
the  geometrical  problem  of  finding  the  cone 
parameters  and  thus  finding  the  sniper  direction. 
The  problem  is  essentially  nonlinear  and  quite 
difficult  if  only  the  shock  wave  arrival  times  are 
known  at  the  given  microphone  positions.  To  make 
easy  the  computation  task,  one  can  propose  here  the 
directional  measurements  realized  by  two 
directional  microphones  placed  in  certain  distance 
from  each  other.  A  number  of  additional 
omnidirectional  measurements  would  improve  the 
localization  accuracy  and  assure  the  solution 
uniqueness.  These  auxiliary  microphones  could  be 
planted  ’in  field’  over  the  safe  -  guarded  area, 
transmitting  their  raw  observations  to  the  computer 
system  by  standard  radio-links. 

It  is  shown  that  only  two  directional  measurements 
are  sufficient  to  evaluate  the  cone  parameters.  The 
angular  accuracy  of  the  measurements  is  expected 
much  lower  than  the  temporal  accuracy  of  the 
shock  wave  arrival  time  measurements,  which  can 
be  assumed  exact.  A  number  of  omnidirectional 
measurements  (yielding  exact  data  of  the  shock 
wave  arrival  times)  help  us  to  overcome  the 
problem  with  inaccurate  directional  measurements 
that  can  contribute  much  to  the  localization 
accuracy. 

2.  The  shock-wave  geometry 

Assuming  constant  velocity  v  and  a  straight  bullet 
path,  the  generated  weak  shock  wave  in  air  has 

a  form  of  a  cone  Sx  (Fig.  1)  of  axis  k  (the 
normalized  vector  directed  against  the  bullet 
velocity,  toward  the  sniper)  and  a  tip  Ox ,  the 


position  of  the  bullet  tip  at  the  given  observation 
time  tx. 


Fig.  1.  The  shock  wave  cone  at  the  observation  time  tx 
when  the  supersonic  bullet  tip  position  is  Ox . 

The  acoustic  signal  of  the  shock  wave  propagating 
with  velocity  C  in  the  direction  normal  to  the  shock 
wave  cone,  arrives  to  the  observers  positioned  at 
r  (in  given  cartesian  coordinates)  at  the  shock 

wave  arrival  times  tt .  All  the  observers  residing  on 
the  shock  wave  cone  Sx  would  detect  the  wave  at 

the  same  time  tx  (simultaneously).  Here  we  only 
consider  the  shock  wave  generated  by  the  bullet  tip 
moving  with  the  supersonic  velocity  V  >  c .  It 
corresponds  to  the  front  edge  of  the  acoustic  N  - 
shaped  signal  [7]. 

At  the  time  t2  >  tx ,  the  shock  wave  cone  tip  moves 
to  the  point  02  : 

02=0i-vk(r2-/i),  (l) 

and  the  cone  broadens  by  the  distance  d  normal  to 
the  cone  S  (Fig.  2), 

d  =  c(t2-tl),  (2) 

where  C  is  the  sound  velocity  in  air  (constant  in  the 
assumed  homogeneous  air).  The  shock  wave 
conical  angle  is 

sini9  =  \c{t2 -fj)]  /  \y{t2-tx )]  =  c/v,  (3) 

assuming  the  supersonic  bullet. 

The  important  conclusion  results  from  the  above 

that  shifting  the  observation  point  r2  (where  the 
shock  wave  arrives  at  the  time  t2 )  by  the  distance 
d  ,  Eq.  (2),  against  the  outward  normal  n2  to  the 
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observed  shock  wave-front  (the  cone  S2),  places 
the  point  r2  on  the  first  cone  Sx  (Fig.  2);  generally 

Yi=Yi-cni{ti-tx),  (4) 

for  zth  microphone  measuring  the  shock  wave 
arrival  time  tt  at  position  r  . 


Fig.  2.  The  shock  wave  expands  over  time  with  the  sound 
velocity  in  air  C  . 

3.  Directional  measurements 


Let  the  two  directional  microphones  placed  at  T\ 
and  r2  detecting  the  shock  wave  arrival  times 
tx ,  t2,  and  simultaneously  the  sound  propagation 

directions  nx  and  n2  (which  are  the  outward 
normals  to  the  shock  wave  cone),  respectively.  It  is 
shown  below  that  these  two  measurements  suffice 

to  evaluate  1)  the  cone  axis  k  and  2)  its  tip  Ox ,  as 

well  as  3)  the  conical  angle  1 9 ,  that  is  4)  the  bullet 
velocity  V ,  and  finally  5)  its  path  in  space 

determined  by  {(9, 

Consider  a  line  described  by  its  point  r  and  the 


vector  ni  along  it.  The  line  is  the  sound  ray 


generated  at  Pt  and  arriving  at  r  at  the  time  tt 
(Fig.  3). 


Fig.  3.  The  cone  axis  k  crosses  two  normals  to  the  cone 
surface  flx  2  at  the  same  angle  i9  . 


The  bullet  path  described  by  the  point  Ox  and  the 
vector  k  crosses  the  rays  ( r  ,  n[ )/  i  =  1;  2,  at 
points  Pt  at  the  same  angle 

3 '  =7t  12-3,  (5) 


as  it  results  from  the  geometry  shown  in  Fig.  2. 
This  yields  the  following  vector  equations  (a  dot 
means  the  scalar  product) 

r-«,n,  =p,,/  =  l,2, 
k  =  ±(P1-P2)/||P1-P2||,  (6) 

k  n1  =  kn2  <  0, 


which  the  last  inequality  helps  us  to  choose  the 
correct  sign  to  k ;  OCi  are  unknown  constants 
(scalars).  On  the  strength  of  Eq.  (2), 

(**2  -  P2  )  *  »2  -  (ri  -  P1  )  •  ni  =d  =  c{t2-tl),(l) 


and  the  second  of  Eqs.  (6)  multiplied  by  ^Px  —  P2 1|  ; 

(r,  -  ax n1  -  r2  +  a2n2 )  •  n1  =  (rx  -  axnx  -  r2  +  ^2112 )  •  n2 

(8) 


one  can  obtain: 


n.  _(r2-ri)-n2+c(^-^2) 
^1  t  ’ 
1  —  n1  •  n2 

l-n1  n2 

sin  i9  =  -n,  •  k  >  0,  /  =  1,  2  , 


(9) 


where  fli ;  k  are  normalized  vectors:  7^.  =  7^. . 
nj  =  1 ,  similarly  k  •  k  =  1 .  The  solution  does  not 
exist  if  nx  •  n2  =  1 ,  which  particular  case  (where 
both  reside  on  the  same  generatrix  of  the  cone)  is 
neglected  in  this  study. 
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The  cone  tip  Ox  can  be  evaluated  from  the  right 
triangle  ( Ox  ,Px,r^),  Fig.  4,  where  both  3  and 

||Pi-ri||  =  (ri-pi)-ni  =  ax  >  (10) 


are  known,  yielding 


Ol=rl+  ax 


IV 


a y 


nl  k  J 

v  =  -c/(nj-k)  . 


=  P,  -  k  ? 
sini9 


(11) 


This  concludes  the  searched  solution. 


Now  consider  an  omni-directional  measurement  of 
the  shock  wave  arrival  time  ^3by  a  microphone 

placed  at  point  r3  (Fig.  4).  It  is  assumed  here  that 

the  same  shock  wave  (the  same  bullet)  is  detected 
by  the  directional  microphones  yielding  the  data 

(r ,  nt ,  t . ),  i  —  1, 2,  and  the  evaluated  shock  wave 
tip  Ox  at  the  time  tx .  The  shock  wave  cone  tip  at 
the  time  1 3  is  03 : 

03  =01-vk(f3-f1),  (13) 


according  to  Eq.  (1).  Naturally,  the  cone  has  the 
same  conical  angle  3  =  arcsin(c  /  v) ,  hence 


Fig.  4.  The  consistency  condition  of  omni-directional 
measurement  at  r3  results  from  comparison  of  the  conical 

angle  3 . 

4.  The  measurement  consistency 
criterion 

Having  three  or  more  directional  measurements  of 
the  same  shock  wave,  the  same  £,vand  other 
parameters  should  result  from  Eqs.  (9),  (11), 
evaluated  from  different  pairs  of  data  but  including 

the  same  (ri,ni,ti)  .  For  example,  Eqs.  (9)  applied 
to  the  pairs  (1;  2)  and  (2;  3),  should  yield  the  same 
a2 

(ri  - r2 ) ’ ni  +  c  (*2-Q_(r2-r3)-ii3+c  (;2-;3) 

1  —  n1  •  n2  1  —  n1  •  n3 

(12) 

in  order  to  obtain  the  same  P2  and  thus  the  same 
k  .  This  consistency  condition  (applied  for  different 
CCi )  will  be  exploited  later  in  order  to  reduce  the 

measurement  errors  causing  the  above  equation  to 
fail. 


(r3-Q3)  k 

Ik-Osll 


=  cos  >9. 


If  measurements  are  exact  then 


(14) 


£  = 


(r3-Q3H 

||r3  °3 1 


+  k  nx 


(15) 


equals  zero  (k  nx  =  —  silli9)  ,  otherwise  t,  ^  0 
indicating  the  incorrect  measurement  data.  The 
condition  ^  —  0  is  the  consistency  condition  of 
omni-directional  measurement  data  with  respect  to 
the  pair  of  directional  measurements  which  can 

include  certain  errors  concerning  the  directions  n  • . 

According  to  the  earlier  assumption  that  the  arrival 
time  is  measured  exactly,  the  omnidirectional 
measurements  are  considered  exact. 

Introducing  matrix  notations  where  vectors 

ri,Oi,k, nt  are  row  matrices  and  r  ,k  etc.  are  their 

transposed  (column)  matrices,  the  above  equation 
can  be  conveniently  rewritten  in  the  form 


^(3)=, 


I z  (l-k'  k )  z ' 


z  z 


+  k  n\,  z  =  r3  - O, 


(16) 

(/is  a  unitary  matrix).  All  the  earlier  vector 
equations  can  be  rewritten  and  evaluated  in 
a  similar  manner,  and  the  same  notations  are 
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applied  to  both  vectors  and  the  corresponding 
matrices. 

The  measurement  consistency  condition  —  0  will 
help  us  reduce  the  directional  measurement 
inaccuracy.  Assuming  r ,  known  exactly  and 

admitting  certain  errors  8i  in  the  measured  shock 
wave  propagation  direction  ili ,  this  condition 
yields  an  equation  for  8i .  Several  omni-directional 

measurements  (at  points  r3  ,  r4  )  are  necessary  to 
obtain  a  sufficient  number  of  equations  in  order  to 
evaluate  all  the  components  of  vectors  8t  of 

interest.  Unfortunately,  Eqs.  (16),  (9)  and  others  are 
highly  nonlinear  and  their  solution  may  not  be 
unique,  in  general.  It  is  assumed  here  that  the 
measurement  can  be  only  slightly  inaccurate  which 
allows  us  to  apply  the  perturbation  analysis  with 

respect  to  8t .  The  resulting  linear  equations  for  the 
measurement  errors  £fcan  be  easily  solved.  It  is 

a  matter  of  numerical  testing  how  large  8t  can  be 
admitted  to  obtain  convergent  solution  for  given 
microphone  positions  r  ,  and  the  bullet  miss- 
distance  and  direction  —  k  . 

5.  Perturbation  analysis 

Assuming  the  measured  (normalized) 

direction  ni  +  8{  instead  of  the  correct  ni ,  it  is 
evident  that  the  equation 

Srnt=  0  (17) 

results  from  the  normalization  condition  ( ni  +  8t )  . 
{ni  +  <^-),  neglecting  higher  order  terms.  Eq.  (17) 
shows  that  8t  has  only  two  independent 

components  orthogonal  to  «. .  We  may  choose 
them  in  directions  of  two  orthogonal  vectors: 


=  d,et  , 

where  dt  is  the  row  matrix  with  two  components 
fully  characterizing  the  measurement  errors  (for  the 
already  chosen  et). 

In  perturbation  analysis  dt  is  infinitesimal,  but  in 
real  computations,  for  8i  small  but  finite,  the 

corrected  vectors  ni  <—  ni  +  St  must  be  always 

normalized  in  order  to  keep  the  earlier  equations, 
like  Eq.  (15),  valid.  The  perturbation  analysis  (for 

infinitesimal  dt )  of  Eqs.  (9)  yields  in  the  matrix 
notations  and  using  the  summation  convention: 


ii 

L 

axn\ 

a2n\+(rl  -r2)' 

axn}2+{r2  -rj 

'  a2n\ 

(20) 


where  ax  2  are  unperturbed  scalars  evaluated  from 
Eq.  (9)  within  zero-order  accuracy  with  respect  to 
8 t .  To  indicate  the  set  of  data:  (^nxdvr\\  and 

(/22,£2,r2),  used  for  evaluation  of  these 
coefficients,  the  superscript  (1;  2)  will  be 

introduced  in  the  subsequent  analysis  like 

ax  ,i  =  1,2  . 

Similarly,  one  can  obtain  the  perturbation  equations 
for 

8 Oi  -  SPji , 8 k  =  Sjkj ,  and  St;  —  Spj 

(21) 

(the  unperturbed  t;  has  zero  value).  Note  that 
again,  Sk  is  an  orthogonal  vector  to  k  : 


eT  =  nt  x  r„ 
e.2)  =  n:  x e,.(1) , 


(18) 


again  normalized  after  evaluation  of  the  vector 
products  denoted  here  by  x.  In  matrix  notations: 


(Sk)k'=  0  (22) 

because  of  normalization  of  k  .  It  can  be  evaluated 
from  the  perturbation  of  Px~  P2,  Eqs.  (6);  the 
perturbation  of  t;  is  evaluated  from  Eqs.  (15) 
provided  that  the  perturbations  to  Ot ,  Eqs.  (11)  - 
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(13)  are  evaluated  at  first.  The  explicit  formula  are 
too  long  to  be  presented  here;  note  only  the 

introduced  perturbation  matrices  O // ,  A+ ,  ^  -  in  Eqs. 
(21)  which  will  be  applied  in  further  analysis. 

6.  The  system  ”2+4” 

Two  directional  measurements  introduce  four 
unknown  errors:  dt  ’  ,Z  =  1,2  (each  §i  has  two 

independent  components  dt  ’  ).  Four  independent 

conditions  from  the  measurement  consistency 
conditions,  for  instance,  are  needed  to  evaluate 

,  /,  j  =  1,2  and  to  retrieve  the  correct  values 

of  ni .  It  shows  that  two  directional  and  four  omni¬ 
directional  (introducing  no  extra  errors) 
measurements  of  the  shock  wave  generated  by  the 
passing  by  bullet  are  sufficient  for  evaluation  of 

correct  n[  in  spite  of  the  measurement  errors.  The 
correct  values  of  n[  allow  one  to  evaluate  the 

correct  bullet  path  parameters,  k ,  Ox  and  V ,  that  is 
needed  to  obtain  correct  sniper  localization.  One 
can  conclude  that  the  directional  measurements 
yield  only  the  first  guess  of  these  parameters  to  ease 
the  corresponding  computation  task  based  on  the 
perturbation  analysis. 

Another  measurement  systems  can  be  proposed  as 
well.  For  example,  three  directional  measurements 

(ni9ti9i  =  1,2,3)  introducing  six  unknowns 
[d^l,2\i  =  1,2,3),  and  three  omni-directional 

ones  (at  different  ri9i  =  4,5,6)  which  can  be 

exploited  for  formulation  of  six  consistency 

(u) 

conditions  A56  ,  evaluated  using  different  pairs 

of  directional  data:  {nx,n^)  and  (^2,^3),  for 
instance.  The  other  possibility  is  to  formulate  three 
consistency  conditions  like  in  Eq.  (12)  for  OCi , 

appended  by  three  chosen  to  obtain  the  best 

conditioned  system  of  equation.  Yet  another  system 
uses  four  directional  measurements  and  two 
omnidirectional  ones;  they  will  be  discussed  below 
in  some  details.  Note  however  that  directional 
measurements  are  much  more  expensive  than  the 
omnidirectional  ones  thus  the  first  above  mentioned 
system,  ”2+4"  is  preferred. 

At  first  glance,  the  system  of  four  directional 
measurements  seems  self-correcting  without 
omnidirectional  measurements.  Namely,  one  can 
formulate  the  sufficient  number  of  consistency 


conditions  like  Eq.  (12)  using  different  pairs  of  the 
directional  measurements  only.  Regretfully,  the 
rank  of  such  system  appears  to  be  only  six 
indicating  that  two  another  equations  are  necessary, 
namely  resulting  from  independent  omnidirectional 
measurements. 

The  system  of  equations  resulting  from  two 
directional  measurement  data:  ,  ni  +  8i  at 

positions  r,  /  =  1, 2,  and  four  omni-directional 
measurements:  t.  at  different  =  3,..., 6, 

results  from  four  consistency  criteria  Eq.  (16). 
Explicitly,  according  Eq.  (21): 

=  xj,i  =  l,2,y  =  3, 4, 5, 6,  (23) 

what  can  be  further  transformed  using  Eq.  (19)  to 
obtain  the  complete  system  of  equations  for  the 

unknown  d  ^ ;  z,  j  =  1,2  : 

=wn =w  <24> 

i,l=\ 

where  the  values  of  X.  are  ^  evaluated  from  Eq. 
(16)  using  the  measurement  data  (that  is  the  values 
of  ni  including  certain  error  8i  what  causes 

This  solved,  yields  the  measurement 
errors  8t  which  subtracted  from  the  measured  data 
yield  the  correct  directions  Tli .  In  practice,  8i  are 
not  infinitesimal  and  the  values  of  8-  '^  are 
evaluated  from  inaccurate  values  ni  —  n[  +  8X . 
Although  the  nt  +  8X  is  considered  closer  to  the 

correct  ni ,  it  is  evident  that  the  correct  solution  can 

be  obtained  repeating  the  calculations  in  a  recursive 
manner.  If  convergent,  they  yield  the  searched 

correct  lfli ,  and  finally  the  correct  bullet  path 
parameters,  particularly  the  most  important  k  . 

7.  Numerical  example 

The  first  considered  system,  "2+4, 1 "  is  the  cheapest 
one.  Numerical  results  show  that  happily,  it 
performs  also  better  than  others,  yielding 
a  convergent  system  of  equations  for  larger 
domains  of  bullet  path  parameters  and  larger 
directional  measurement  inaccuracy. 
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In  the  numerical  example  presented  here  for  the 
system  "2+4",  the  microphones  are  placed  on  the 
ground  (this  is  also  the  sniper’s  post  level)  and 
distributed  over  the  protected  area  about  20m  long 
(Fig.  5;  squares  represent  directional  microphones, 
circles  -  omnidirectional  ones).  The  bullet 
missdistance  is  assumed  2m  above  the  ground,  and 
its  velocity  is  3 C  .  The  directional  inaccuracy  is 
modeled  by  performing  calculations  for  20  random 
directions  ili  +  8X  within  1  %  limit  off  the  correct 


Fig.  5.  The  simulation  of  the  system  ”2+4"  for  1% 
measurement  inaccuracies  of  directional  measurements. 
Axes  units  are  l[m].  Dot  lines  shows  successful 

evaluation  of  the  sniper  direction  k  ,  and  solid  lines 
indicate  cases  of  not  convergent  iterations. 

If  all  calculations  converge  to  correct  ni  than  the 

corresponding  bullet  path  is  plotted  with  dash, 
otherwise  with  solid  line.  The  example  shows  that 
the  system  fails  in  only  few  cases  of  sniper’s  fire 
aimed  at  different  points  from  different  sniper’s 
posts. 

More  extensive  simulations,  beyond  the  scope  of 
this  paper,  for  realistic  cases  of  the  microphone 
distribution,  sniper’s  positions  and  bullet  path 
orientation  with  respect  to  the  protected  area,  would 
reveal  the  true  value  of  the  above  proposed 
systems. 

Having  the  bullet  path  evaluated  (characterized  by 
Ox,k,V ),  one  can  easily  exploit  the  other  acoustic 
information  about  the  fire  -  the  acoustic  signal  of 
the  muzzle  blast.  Assuming  a  small  miss-distance 
of  the  sniper’s  fire  and  sufficiently  large  distance 
( L )  to  the  sniper,  the  muzzle  blast  signal 
propagates  nearly  along  the  bullet  path.  This  gives 
the  approximation  concerning  the  time  difference 

between  the  measured  shock  wave  (tx)  and  muzzle 
blast  ( t0  )  arrival  times: 


L/c-L/v  =  t0-tx  (25) 

(tx,t2,  ...  are  assumed  close),  from  which  one 

easily  find  L  and  thus  the  sniper’s  position 
measured  along  the  evaluated  bullet  path. 

8.  Conclusions 

Concluding,  the  acoustic  system  is  proposed  for  the 
localization  of  the  sniper  position  fully  exploiting 
the  most  reliable  [8]  and  impossible  to  mask 
information  delivered  by  a  supersonic  bullet. 
Making  two  or  more  directional,  and  a  number  of 
supplemental  omni-directional  measurements  of  the 
shock  wave  signal  with  adequate  accuracy  (note 
that  this  is  a  ’single-event’  measurement  that  cannot 
be  repeated  to  improve  it),  one  can  evaluate  the 
bullet  path  parameters  (particularly  the  most 
important  bullet  velocity  direction  -  k  pointing  at 
sniper’s  post)  and  the  bullet  velocity  with  improved 
accuracy.  The  proposed  system  has  the  ability  to 
correct  the  directional  measurement  inaccuracies. 
The  numerical  examples  based  on  perturbation 
linearization  of  highly  nonlinear  equations 
governing  the  geometry  of  the  considered  problem 
show  that  this  self-correcting  ability  works  well  for 
the  measurement  errors  as  large  as  1%,  which  can 
be  expected  technically  feasible.  Fully  nonlinear 
analysis  would  certainly  even  lower  this 
requirement  admitting  still  larger  directional 
measurement  inaccuracy  and  making  the  presented 
concept  of  the  sniper  localization  even  more 
attractive  for  implementation,  saving  precious  life 
of  peace-keepers  and  innocent  civilians  from  the 
snipers’  thread. 

Instantaneous  evaluation  of  the  sniper’s  direction 
(-  k )  is  essential  for  this  reason  also  that  it  may 
enable  us  to  apply  the  other  countermeasures.  For 
example,  properly  directed  optical  (infrared)  sensor 
can  pick-up  the  cloud  of  hot  muzzle  gases  detecting 
precisely  the  sniper’s  post.  In  any  cases,  fast 
response  is  necessary  to  prevent  the  change  of  the 
sniper’s  post;  perhaps  automatic  counter-fire  is 
necessary.  Having  known  of  such  countermeasures, 
the  sniper  would  try  to  act  fast,  which  would  surely 
contribute  to  degradation  of  his  fire,  making  it  less 
lethal  in  any  case. 
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